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ABSTRACT 


In our solar system, Mars-sized protoplanets frequently collided with each 
other during the last stage of terrestrial planet formation called the giant im¬ 
pact stage. Giant impacts eject a large amount of material from the colliding 
protoplanets into the terrestrial planet region, which may form debris disks with 
observable infrared excesses. Indeed, tens of warm debris disks around young 
solar-type stars have been observed. Here, we quantitatively estimate the to¬ 
tal mass of ejected materials during the giant impact stages. We found that 
~0.4 times the Earth’s mass is ejected in total throughout the giant impact 
stage. Ejected materials are ground down by collisional cascade until micron¬ 
sized grains are blown out by radiation pressure. The depiction timescale of 
these ejected materials is determined primarily by the mass of the largest body 
among them. We conducted high-resolution simulations of giant impacts to ac¬ 
curately obtain the mass of the largest ejected body. We then calculated the 
evolution of the debris disks produced by a series of giant impacts and depleted 
by collisional cascades to obtain the infrared excess evolution of the debris disks. 
We found that the infrared excess is almost always higher than the stellar infrared 
flux throughout the giant impact stage (~100 Myr) and is sometimes ~10 times 
higher immediately after a giant impact. Therefore, giant impact stages would 
explain the infrared excess from most observed warm debris disks. The observed 
fraction of stars with warm debris disks indicates that the formation probability 
of our solar system-like terrestrial planets is approximately 10%. 

Subject headings: planets and satellites: formation - planets and satellites: terrestrial 
planets — protoplanetary disks 
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Introduction 


Collisions among Mars-sized protoplanets have been though t to be co mmon events 


during terrestrial 


Agnor et al. 


planet formation in our solar system (e.g., 


Chambers fe Wetherill 


1998 


19991) . and these collisions are called giant impacts. In the solar system, several 


Mars-sized protoplanets formed in the terrestrial plan et reg ion 


through the successive accretion of planetesimals (e .g., 


Wetherill 


1AU from the Sun) 


1985 


Kokubo fe Ida 


1998 


Kobavashi et al. 

2010; 

Kobavashi & Dauohas 

2013) 


20131) . Giant impacts among protoplanets 


are triggered by the depletion of the gas in the protoplanetary disk, which leads to the 
instability of the orbits of the protoplanets. The stage in which many giant impacts occur is 


Known as me g 

Jacobson et al. 

lanr n 

2014; 

npaer srage ana iasxs a 

Chambers & Wetherill 

)proxi 

1998) 


1998 ). This is consistent with 


the formation of the Moon, as inferred from its geochemistry flTouboul et ah 
20081). 


he timing of 


20071; 


Hallidav 


Even in extrasolar planetary systems, giant impact events are thought to be common. 
For example, short-period super-Earths have been thought to grow through collisions among 
Mars- and Earth-siz ed proto planets that migrate inward to the vicinity of the central stars 


(Terquern & Pa palo izou 


2007 b These giant impacts in extrasolar planetary systems after 


the depletion of the gas in the protoplanetary disk can produce the observed architecture 


of super-Earth system s where short-period super- 


(jOgi hara fe Ida 


2009 ). 


Ogihara et al. 


([2013 


larths are not in mean-motion resonance 


2 0141 ) also investigated the possibility that 


super-Earths grow through giant impacts outside a hot Jupiter. They found that growing 
super-Earths push the hot Jupiter inward, and eventually the hot Jupiter is swallowed into 
the central star, which may explain the lack of super-Earths in systems with a hot Jupiter. 

Giant impacts significantly influence various features of planet ary sy stems, such as 


the number, mass, and spin of terrestrial planets (e.g., 


Agnor et al 


1999 


Kokubo et al. 
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20061). Gian 
Moon (e.g., 


impacts are highly energetic even ts and are responsible for the creation of t he 


Canup 


2004 : 


Genda fe Ab 


2003afl . the martian s atellites ( Ci 


and planets with extremely large cores such as Mercury (e.g., 


Benz et al 


non et al. 


2015), 


2007). Moreover, 


giant impact s are closely related to t he thermal states of terrestrial planets, such as magma 


oceans (e.g., 


Tonks fc Melosh 


planet atmospheres flGenda fe Ab 


19921 


Harnano et al. 


2003b, 


2005 ). 


201311 . and the origins of terrestrial 


According to recent simulations of giant impacts, there are several types of collision 


outcomes, including 


(lAgnor & Asphaug 


2004 


i — - r>—n? — 

Leinhardt & Stewart 

1 ? r ^ 

2012; 

Genda et al. 

2012; 


isruptive collisions 


Sekine fe Genda 


2012). Even in the case of a low-velocity giant impact (vi mp ~ v esc , where Ui mp and v e 


are the impact velocity and the two-body surface escape velocity, respectively), such as 


the Moon-forming giant impact, approximate 


escapes from the Earth’s gravity (e.g., 


Canup 


y 1% of the Earth’s mass is ejected and 


20041 1. Moreover, if Mercury’s large core is 


the result of gi ant impact(s), a s i gnificant amount of mantle ma terial is ejected from the 


proto-Mercury ( Benz et al. 


1988 


2007; 


Asphaug & Reufer 


2014). Therefore, it is expected 


that a large amou nt of material is ejected into the terrestrial plane 


impact stage (e.g., 


Leinhardt fe Stewart 


2012 


Stewart fc Leinhardt 


region during the giant 


201211 . In this paper, we 


call the material ejected by a giant impact “fragments” regardless of their sizes and masses. 

Recently, thanks to infrared space telescopes such as Spitzer and AKARI , tens of warm 
debri s disks around solar-type (FGK) stars with a ges of 10 7 — 10 8 years have been reported 


(e-g., 


Zuckerman et ahl 


2011 


Fujiwara et al. 


2013.1). The temperatures of these debris disks 


have been estimated to exceed 150 K from the spectral energy distribution including 
infrared excess. These warm debris disks are estimated to be located roughly 1 AU to 
several AU from the central stars, which corresponds to the terrestrial planet region in the 
solar system and the super-Earth-forming region around some extrasolar stars. Based on 
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their stellar ages and locations of the debris disks, the relation b etwee n these warm debris 


Jackson & Wvatt 

to 

o 

to 

Mclis et al. 

2012) 


Weinberger et al. 


2011 


Jackson .& W y att 


(2012) considered the Moon-forming giant impact. This impact is 


an oblique collision between Earth- and Mars-sized protoplanets with a low impact velocity 
(~ n eS c), which led to the formation of a proto-lnnar disk, from which the Moon was 
made, around the Earth-sized protoplanet. They focused on the fragments produced by 
the Moon-forming giant impact, and performed numerical simulations of the orbital and 
collisional evolution of the ejected fragments around the Sun. They adopted 1.6% of M e 
(the mass of the Earth) as the total mass of the eje cted fragments , w hich was taken from 


the numerical result of the giant impact simulation flMarcus et ah 


2009l h and 0.2% of M e 


as the largest fragment, which was taken from the largest clump in the proto-lunar disk 
produced by the Moon-forming giant impact (Canup 2008). The lifetime of a debris disk 
strongly depends on the mass of the largest fragment. Although the mass of the largest 
clump in the proto-lunar disk is not exactly the same as that of the largest fragment that 
is ejected from the potential well of the Earth-sized protoplanet, they estimated that the 
debris disk produced by the single Moon-forming giant impact has a lifetime of ~ 30 Myr. 
They concluded that this debris disk can be observed by infrared excess at 24 /im. 


Many giant impacts occur during the giant impact stage, and there are several types 
of collision outcomes other than Moon formation. To determine the infrared excess from a 
debris disk throughout the giant impact stage, all types of giant impacts that take place 
in the giant impact stage should be taken into account. Moreover, because the total mass 
of the fragments and the mass of the largest fragment ejected by each giant impact are 
the important parameters in determining the lifetime of the debris disk, these parameters 
should be determined via high-resolution giant impact simulations. 
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In this paper, we calculate the infrared excess from a debris disk produced by multiple 
giant impacts in the giant impact stage. For this purpose, we need to determine the 
number and timing of giant impacts, the outcomes of each giant impact, and the collisional 
evolution of the fragments. The giant impact events in the giant impact stage occur 
because of the chaotic orbital instability of protoplanets, which must be investigated by 
IV-body simulations taking the mutual gravity int o acco unt. Using the results of iV-body 


simulations of the giant impact stage by 


Kokubo & Gendal (120101) and the smoothed particle 


hydrodynamics (SPH) simulations of giant impacts by Genda et al. (2012), we calculate the 
total ejected fragment mass during the giant impact stage, as described in Section [21 The 
fragments collide with each other, which leads to produce smaller fragments. This collisional 
cascade grinds down fragments to micron-sized grains, which are blown out by radiation 


Kobavashi & Tanaka 

2010 

). However, the SPH simulations by 

Genda et al. 

(2012) 


not performed at a sufficiently high resolution to accurately obtain the mass of the largest 
ejected fragment. In Section EJ we newly conduct high-resolution SPH simulations of giant 
impacts. In Section [I] we investigate the collisional evolution of the fragments using the 
results obtained in Section El and we calculate the mass evolution of debris disks. Using 
these results, we show the infrared flux evolution during the giant impact stage in Section 
El Finally, we discuss our findings in Section El 


2. Total Mass Ejected by Giant Impacts during Giant Impact Stage 

In this section, we quantitatively estimate the total mass of fragments ejected by 
giant impacts tha t occu r duri ng the giant i mpac t stage by using the results obtained by 


Kokubo & Genda 

(2010 

) and 

Genda et al. 

(2012) 





















Kokubo & Gendal (2010) performed A-body orbital calculations of protoplanets in 
the giant impact stage np to 200 Myr, taking into account both hit-and-run and merging 
collisions. For the initial conditions, 16 Mars-sized protoplanets with a total mass of 2.3 M ( 
are located from 0.5 to 1.5AU. The orbital separations of protoplanets are given by_10 


mutu al Hill radii in accordance with the results of planetesimal accretion ( Kokubo fe Ida 


1998 ). The initial eccentricities and inclinations of protoplanets are s et to h ave a R ayleigh 


istribution with dispersions of 0.01 and 0.005 rad, respectively. In 


Kokubo fe Genda 


(201Q), 50 runs with different initial angular distributions of protoplanets were performed. 


As the results of statistics, 1211 giant impacts occur in 50 runs. This corresponds to an 
average of 24 giant impacts per run. The number of finally formed terrestrial planets is 
3.6 ± 0.8, and the final giant impact occurs at 73 ± 74 Myr. These numerical simulations 
provide us the time, location, and impact parameters (the masses of the two colliding 
protoplanets, impact velocity, and impact angle) of all giant impacts that take place during 
the 50 runs of the giant impact stage. 


Genda et al. 


(20121) performed SPH simulations of giant impacts between protoplanets 


with a 30wt% iron core and a 70wt% silicate mantle. More than 1000 impact simulations 
for various impact parameters (the mass ratio of the colliding protoplanets, impact velocity, 
and impact angle) were systematically car ried out, and the merging criteria of protoplanets 


were investigated. Although 


Genda et al 


(2012) did not focus on the ejected fragments, 


in thi s paper , we calculate the total mass of the fragments based on the data obtained by 


Genda et al. 


(2012). Figured] shows an example of the total mass of the fragments ejected by 


a single giant impact between protoplanets of the same size with mass 0.1M® as a function 
of impact velocity, v imp , and angle, 6. Collisions with lower angles (near-head-on collisions) 


and higher impact velocities tend to produce a larger amount of fragments. We only count 
the ejected fragments that are not gravitationally bounded to the colliding protoplanet(s), 
and thus the fragments orbiting around the protoplanets, such as a proto-lunar disk, are 





















not included in this mass. We tabulate the data of the ejected mass M e j e for various impact 
parameters. 


Using the impact parameters obtained by Kokubo and Genda (2010) and the table 
data of M e j e made here from the results obtained by Genda et al. (2012), we can calculate 
M e je for all giant impacts that take place in the 50 runs of the giant impact stage. Figure 
[2] shows the total mass of the ejected fragments, MU), in each run. The average total 
mass over the 50 runs is MU) = 0.42M e , which corresponds to 18% of the total mass 


of the system (2.3 


Stewart fe Leinhardt 


and is consistent with the previously estimated value of 15% by 


(2012). The maximum and minimum total masses are MU) = 1.76M a 


and 0.09M ffi and occur during Runs 21 and 26, respectively. On average, each giant impact 
produces fragments with a total mass of 0.02M e . Although the ejected mass by a single 
giant impact is much smaller than the mass of the terrestrial pla nets, MT is comparable to 
or greater than the typical mass of warm debris disks (e.g., 


Fujiwara et al. 


2012). ft should 


be noted that the ejected mass is reduced by collisional evolution in the debris disk. We 
investigate the evolution of the debris disk in Section |4j 


The values of MU) estimated here is somewhat overestimated, because 


Kokubo fc Genda 


(120101 ) does not take into account the mass loss by giant impacts in their IV-body calculations. 


We used the masses of perfectly merged protoplanets to estimate each M e j e . Especially, for 
the cases with high MU) such as Runs 8 and 21, the protoplanets should be significantly 
smaller than perfectly merged protoplanets in the latter part of the giant impact stage. We 
will discuss this issue in the next section. 
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3. Giant Impact Simulations 


The fragments ejected by a giant impact evolve via successive collisions among them. 
Small fragments are ground down to the blow-out size, which is ~ 1 fi m around solar-type 
stars, on a short timescale, whereas large fragments reach the blow-out size later via 


collisional cascade. Therefore, the decay timescale of the fragments is dete rmined primari 


by the mass of the largest fragment and the total ejected mass (IKobavashi fe Tanaka 


2010 ). 


y 


In this paper, the largest fragment does not refer to th e post- impact protoplanet(s), but 


the largest body ejected by a giant impact. Although 


Genda et al 


(2012) performed giant 


impact simulations using 20,000 SPH particles, their numerical resolution is insufficient to 
determine the mass of the largest fragment. Moreover, it is difficult to ma ke a table of 
the mass of the largest fragment using the data obtained bv iGenda et al] (2012), because 


the formation process of the large fragments is rather stochastic, as shown below. In this 
section, we perform higher-resolution calculations (100,000 SPH particles) of giant impacts 
that took place in Run 1, 21 , and 26 from the 50 runs, using the exact impact parameters 


of giant impacts obtained by 


Kokubo &; Gendal (120101) to precisely calculate the mass of the 


largest fragment ejected by a giant impact. 

To accurately calculate the mass of large fragments, we use the calculation method 


developed by 


Genda et al. 


(120151) . In this method, we use a friends-of-friends algorithm 


to roughly identify clumps of SPH particles. Then, we iteratively determine if each SPH 
particle is gravitationally bound to these clumps. Finally we determine if these clumps are 
gravitationally bound to each other. In this analysis, we set the minimum number of SPH 
particles for a clump to be 10. If any clumps except for protoplanets are not detected, we 
set ilfirg to be the mass of just one SPH particle. 


Our numerical c ode is the standard SPH method (e.g., 


same as that used in 


Genda et al. 


Monaghan 


1992), which is the 


( 2012 ). We performed the SPH simulations of giant 
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impacts between protoplanets with a 30wt% iron core and a 70wt% silicate mantle. We 


applied the Tillotson equation of state (Tillotson 
of the protoplanets were set to zero. 


1962) in our simulations. The initial spins 


Figure [3] sh ows s napsho 


Run 1 of 


;s o f the numerical simulation for the 11th giant impact in 


Kokubo & Genda (|2010|). After the first contact between two protoplanets, they 


escape from each other (hit-and-run collision). Numerous SPH particles are ejected by this 
collision. In this simulation, fragments with a total mass of 0.05AT® are ejected. Many 
clumps are gravitationally formed from ejected SPH particles. The mass of the largest 
clump (here, we call the largest fragment) is approximately half of the lunar mass, which 
includes approximately 2,500 SPH particles in this simulation. 


We performed numerical simulations of all giant impacts that took place in Runs 1, 21, 
and 26 from the 50 runs, which correspond to the cases of the average (MT) = 0.47M®), 
maximum (A/fW = 1.76M®), and minimum (M*?* = 0.09M®) total ejected masses estimated 
in the previous section, respectively (see Figure [2]). In Runs 1, 21, and 26 of the simulation, 
35, 47, and 17 giant impacts occurred, respectively. Figure [I] shows a plot of the total ejected 
mass M e j e and the mass of the largest fragment Mi rg for each giant impact as a function 
of time. Calculated MA* from these high-resolution giant impact simulations is 0.45A7®, 
1.59AT®, and 0.083A7® for Runs 1, 21, and 26, respectively, which are consistent with 
previously estimated AT*?* from the table data of low-resolution giant impact simulations. 


As we mentioned in the previous section, the values of MgP estimated here is somewhat 


overestimated. To roughly estimate the effect of the mass 


oss from protoplanets, we 


Ko kub o & Genda ( 2010 ) and 


reanalyze A-fT). We use the impact parameters obtained by 
collision outcomes obtained from high-resolution impact simulations, but we take into 
account the effect of the decrease in the mass of protoplanets by giant impacts. The 
decrease in M e j e is also considered according to the mass of colliding protopalnets whose 
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masses are equal to or smaller than those of perfectly merged protoplanets. We find that 
the effect of the mass loss reduces from 0.45M® to 0.38M e for Run 1, from 1.59M® to 
0.82M e for Run 21, and from 0.083M e to 0.081M® for Run 26. The effect of the mass loss 
is significant for Run 21, and is overestimated by a factor of 2. However, significant 
amount of fragments are still ejected during the giant impact stage. On the other hand, 
this effect is small for Run 1 and negligible for Run 26. 

To compare our results with the stellar age in Section [6l we assume that the giant 
impact stage starts at 10 Myr, because the formation of protoplanets followed by formation 
of planetesimals in a protoplanetary disk would take several Myrs, and the giant impacts 


are triggered by the gas depletion of 

protoplanetary disks, whic 

i also takes several Myrs 

inferred from the disk observations ( 

Beckwith & Sargent 

1996; 

Wvatt et al. 

2003 

). As 


shown in Figure HI most giant impacts occur within 100 Myr. In Run 21, a few disruptive 
giant impacts (M e j e > 0.1M ffi ) occur. The mass of the largest fragment M\ vg strongly 
depends on each giant impact event. As shown in Figure [3j clumping is caused by the 


mutual gravity of ejected SPH particles. This process is rather stochastic, and it is difficult 
to estimate M\ rg without numerical simulations. As shown in Figure HI some giant impacts 
result in producing large M\ Ig , and some result in very small M\ vg . We can safely say that 
Mi rg is less than 10% of M e j C in most cases. In Run 26, almost all giant impacts are gentle 
events, and after the last giant impact at 20 Myr, the planetary system is stable. 


4. Lifetime of Debris Disks 

Fragments ejected by a giant impact are further ground down by mutual collisions 
among them until radiation pressure blows away micron-sized or smaller particles. This 
collisional cascade reduces the total mass of the fragments M tot . The mass evolution of 
the fragments ejected by a single giant impact via collisional cascade is given by (e.g., 
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Kobavashi & Tanaka 2010) 

MUt) = 


M l - 


(t > 4n)> 


( 1 ) 


1 + (^ - 4u)/4 P 

where fde P is the mass depletion timescale due to collisional cascade and £ co u is the time of 
the giant impact. The superscript i represents the z-th giant impact. The size distribution 
in the steady-state collisional cascade is analytically derived as Eq. fl5]) in Section [5l because 


Kobavashi & Tanaka 

2010 

). The mass 

ux using the steady-state size distribution around 

M lrg gives t dep as ( 

Kobavashi & Tanaka 

2010 

), 


f dep = 4.2 x 10 6 


M, 


lrg 


V' 64 ( a V 18 (&a/a\ ( e y™ ( M eje y 1 

) ( iau/ (.yyjlyj Uxio^ J years 


6 x 10 22 kg 


( 2 ) 


where a is the semimajor axis at which a giant impact takes place, A a is the range of 
the semimajor axes of fragments, and e is the typical eccentricity of the fragments. This 
formulation of the mass depletion timescale is derived by the analysis including erosive 
collisions caused by low impact energy, which are more effective in the mas s ev olution due to 


the collisional cascade, and assuming the collisional strength derived from 


Benz & Asphau 


( 19991 ) for impacts of rocky bodies at a speed of 5 krn/s. 


The values of M e j e and M\ rg are derived from the SPH im pact s imulations in Section 


[3j and a is given by the location of giant impacts obtained by 


Kokubo & Genda (201 


Because fragments are stirred by the post-collision protoplanet (s), e and Aa are estimated 
from the mass of the protoplanet. We set e = u esc /uK and A a = 10rn, where v esc is the 
surface escape velocity of the protoplanet (r? e sc = \/2GM p / R p , where M p and R p are 
the mass and radius of the protoplanet, respectively), vk is the Kepler velocity of the 
protoplanet y = \/GM s /a, where M s is the mass of the central star; here, we use the 
solar mass), and r H is the Hill radius of the protoplanet (r H = (M p /3M s ) 1 / 3 a). If the giant 
impact is a hit-and-run collision, we use the mass of the larger post-collision protoplanet to 
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estimate e and A a. 


It should be noted that the stirring by a protoplanet increases the random velocities 
of fragments up to v esc because fragments with random velocities larger than v esc effectively 
experience collisions with the protoplanet rather than the stirring. The random velocities, 
comparable to n esc , are much smaller than the Keplerian velocity around 1 AU. Therefore, 
the scattering ejection of fragments from the central star by protoplanets is negligible. 

On the other hand, collisions between a protoplanet and fragments reduce the total mass 
of fragments. The collisional timescale with a protoplanet whose mass is M pro is given 
by (risPcoiQ )" 1 , where n s is the surface number density of the protoplanet, P co \ is the 
collisional probability, and 0 is the orbital frequency. Since the number of protoplanets in 
a fragment annulus is one, n s ~ l/27raAa. The collisional probability is approximated to 

ative velocity between fragments and the protoplanet is v esc 


be 7 tR 2 assuming that t 

re re 

(Greenzweig & Lissauer 

1992) 


onto a protoplanet is given by 


t 

U £LCC ^ 


1 ' 2xl07 fe)" /3 (m) 7/2 ( Ao/ ° 


0.1 


years. 


(3) 


Because t acc is longer than t d e P as shown in Eqs. ([2]) and ([3]), the accretion onto protoplanets 
is negligible. Once the total mass of fragments is smaller than 10 _2 M® due to collisional 
cascade, the accretion onto a protoplanet mainly reduces fragments. However, the disks 
with such small fragment masses are much fainter than those with which we are concerned. 
Therefore, we only consider the depletion of fragments by collisional cascade. 


A collision between protoplanets ejects fragments, and the total mass of the fragments 
decreases by collisional cascade. Another subsequent impact increases the total mass of 
fragments. Because each f dep is generally shorter than the time intervals of the giant 
impacts, for simplicity, we deal independently with the mass evolution of fragments 
produced by each giant impact. The total mass of the debris disk M disk as a function of 
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time throughout the giant impact stage is given by the summation of M t * ot (t) over all giant 
impact events, 

M di!k (i) = £+,,(«)■ (4) 

i 

This independent treatment makes fdep slightly overestimated, because overlapped region 
of debris disks produced by giant impacts has larger M e j e . Figure ED shows the evolution of 
Mriisk for three runs. Here, we assume that the giant impact stage begins at 10 Myr. Each 
spiky structure in the figure corresponds to a giant impact event. The collisional cascade 
following a giant impact reduces M disk , whereas subsequent giant impacts increase M disk . 
As a result, repetitive giant impacts ensure M disk remains larger than 0.01M e for Runs 1 
and 26, and larger than 0.003 M e for Run 21 until ~ 100 Myr. 


5. Infrared Excess of Debris Disks During Giant Impact Stage 


In collisional cascade, the size distribution of the fragments is determined by collisional 
equilibrium. The differential surface number density of fragments with mass m , defined 


as n ft (m)dm, is proportional to m q dm, and the index q is given by (iKobavash i & Tana ka 
2010 ) 

11 + 3 P 

q = - ~ , — , (5) 


6 + 3 p ’ 

where p = d[ln (Q^/u 2 )]/d[lnm], is the specific impact energy needed for the ejection 
of half mass of the colliders, and v is the impact velocity between bodies. More massive 
bodies tend to have larger Q values in the gravity regime (for kilometer-sized or larger 
bodies) because of the strong gravity, whereas Q is constant or has a negative dependence 
on m in the strength regime (for kilometer-sized or smaller bodies). Because v is stirred 
by protop lan ets, v is independent of m. Therefore, using the mass dependence of Q 


(jBen z fe As phau g 


19991 ). we set q = 1.68 for s > s 0 and q = 1.89 for s < s 0 , where s is the 


radius of the body and Sq — 1 km is the approximate value determined from the simulation 
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by Benz and Asphaug (1999). In addition, the micron-sized or smaller grains in the debris 
disks around the solar-type stars are blown out by radiation pressure on a short timescale 
comparable to the Keplerian period (~ 1 year), and hence we set the smallest size of this 
size distribution as 1 /irn. Using this size distribution, we obtain the vertical optical depth 
of the debris disks r from Mauk- 

For the warm debris disks considered in this study, the peak wavelength of thermal 
emission is not much longer than the smallest dust in the debris disks, which is determined 
by radiation pressure blow-out. We thus approximate the thermal emission of dust as 
blackbody radiation. In addition, debris disks are optically thin, because r just after a giant 
impact is typically < 0.1 at the orbit where a giant impact takes place, and r decreases 
quickly with time. Thus, we can express the emission flux with the wavelength A from the 
debris disk as 

Ff k (A) = dj / 2irrr(r)B„(T)dr, (6) 

where D is the distance of the star from the Earth, r is the radial distance of the debris 
disk from the central star, and B u is the Planck function for a given temperature T. The 
temperature is defined as T = 280(r/lAU) 1 / 2 K. Equation (|6]) is integrated over the debris 
disk, typically from 0.3 AU to 2 AU. 

The summation of the disk flux F^ lsk (A) and the stellar flux F* (A) gives the 
observational flux F° bs (A). Figure [6] shows the flux ratio F° bs (24/xm)/F*(24^m) for Runs 
1, 21 and 26, where we use the solar value of F*(24/im). A giant impact increases the flux 
ratio, which is decreased by the subsequent collisional cascade. Multiple impacts lead to 
the spiky evolution of the flux ratio, which is similar to the mass evolution of the debris 
disk shown in Figure [5] The calculated observational flux is comparable to the stellar flux 
until 100 Myr. Our results indicate that bright warm debris disks are generally formed if 
planetary systems have a giant impact stage, as our solar system had. 
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6. Discussion and Conclusions 


In the late stage of terrestrial planet formation in our solar system, Mars-sized 
protoplanets that formed through the accretion of planetesimals collide with each other 
because of the chaotic orbital instability after gas depletion, resulting in the formation of 
Earth and Venus. This scenario is consistent with the core formation timescale of Earth 
(~J.OO Myr) a nd Mars (< 


2008 


Dauphas fe Pourmand 


5 Myr ) obta ined from Hf/W isotope analysis (e.g., 


2011 


Kobavashi & Dauphas 


Hallidav 


20131 ). The f inal giant imp act in 


the formation of Earth may result in the formation of the Moon (e.g., 


Canup 


20041 ). We 


estimated the total mass of the fr agme nts ejected by 


results of 


Kokubo & Genda ( 2010 ) and 


Genda et a! 


giant impacts using the numerical 


(2012). We found that a significant 


amount of fragments (~ 0.4M® in total) is ejected into the terrestrial planet region 
throughout the giant impact stage (~ 100 Myr). Moreover, we accurately calculated the 
mass of the largest fragment using high-resolution SPH simulations of giant impacts. Using 
these results, we obtained the evolution of the debris disks throughout the giant impact 
stage. A single giant impact produces some amount of fragments, and their total mass is 
decreased by collisional cascade among these fragments. A series of giant impacts forms 
debris disks, and the disk’s infrared flux has a spiky time evolution from 10 Myr to at least 
100 Myr (see Figure [6]). 


Our simulation results can be compared with observational results. The observed 24 
/irn flux ratios for warm debris disks (> 150 K) are also plotted in Figure El The debris 
disks caused by giant impacts can explain most of the observed warm debris disks around 
stars with ages of 10 - 100 Myr. However, HD 145263 and HD 113766, which have ages 
of ~ 10 Myr, have extremely bright disks with flux ratios of ~ 100, and our calculated 
evolution of infrared excesses during the giant impact stage cannot readily explain their 
existence. There is a possibility that disruptive giant impacts that occurred at 40 Myr 
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and 150 Myr in Run 21 took place just at the ages of HD 145263 and HD 113766. In 
addition, if vapor condensates produced by high-velocity collisions between protoplanets are 
considered, extremely sharp spikes in infrared excess would be expected. This is because 


vapor materials re-condense into smal 


few cm (e.g., 


Melosh fc Vickerv 


1991 


drop l ets w ith char acter istic sizes of a few mm to a 


Johnson & Melosh 


2012 


2014J ). and these mm-cm 


size spherules grind away very quickly and produce huge infrared excess. However, because 
an excess duration of flux ratios of 100 is very short (< 1000 yrs), the observed probability 
should be extremely low. On the other hand, HD 109085, which is ~ 1 Gyr in age, has a 
debris disk with a flux ratio of ~ 2, which also cannot be explained by the giant impact 
stage, because the final impact usually occurred earlier than 300 Myr in the formation of 
our terrestrial planets. From a theoretical point of view, it is difficult to prolong the giant 
impact stage up to 1 Gyr in our solar system. Another mechanism may be required to 


escribe the c 


([Gomes et ah 


ebris disk formation after terrestrial planet formation, such as the Nice model 


20051 ). in which gas giant planets move after planet formation. 


Debris disks are also formed during the growth of protoplanets through the accretion 
of planetesima 
planetesimals 


Kenvon & Bromley 

2002; 

Kennedv & Wvatt 

2010; 

Kobavashi & Lohne 


2 0141 ). In systems where the giant impact stage is expected, the collisional destruction of 


planetesimals during protoplanet growth forms cold debris disks around stars with ages 
of ~ 10 - 100 Myr. We can distinguish between planetesimal-induced cold and giant 
impact-induced warm debris disks. However, if planetesimal disks are less massive, the 
formation and growth of protoplanets in the terrestrial planet formation region occur 
even at 1 0 - 1 00 Myr, resulting in warm debris disks around stars with ages of 10 - 100 


Myr flKobavashi fe Lohne 


2 014 ). The time evolution of disk fluxes originating from the 


formation and growth of protoplanets is smooth, which is quite different from the debris 
disks caused by giant impacts. We will be able to determine major origin of warm debris 
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disks when the number of sampled debris disks is sufficiently high. Spectral feature of 
warm debris disks produced by giant impa cts wo uld be also different from that produced 


by planetesimal collisions (e.g., 


Lisse et ah 


2009 ). This is because protoplanets should be 


differentiated, and fragments ejected by giant impacts would be rich in crustal component, 
which is different from the composition of undifferentiated planetesimals. 

Our results show that debris disks produced by giant impacts can be observable at 24 
/irn infrared excess throughout the giant impact stage (~ 100 Myr). However, not all stars 
with ages of 10 - 100 Myr have infrared excess. This observational result might suggest 
that the giant impact stage that we consider in this paper is not very common in extrasolar 
systems. Observations of FGK stars with ages of 10 - 100 Myr by the Spitzer space 


telescope have shown t 

hat the fraction of stars with 24 /irn infrared excess 

ranges from 8% 

(Beichman et al. 

2006; 

Carpenter et al. 

2009; 

Trilling; et al. 

2008 

) to 48% ( 

Zuckerman et al. 


2011). Observation of these stars with multiple wavelengths gives the temperature of the 


debris disk. The fraction of stars with the component of a warm debrisdisk (> 150 


among the stars wit 

h 24 / 

36% ( 

Morales et al. 

2011) 


i 24 /irn infrared excess ranges from 13% ( Jackson & Wyatt 2012) to 


<) 


2 0111) . Therefore, the fraction of stars with a warm debris disk is 1 


17% in total, which is also cons istent with the recent observations at 12 /irn by WISE 


space telescope (Kennedy fc Wvatt 


2013|). As described a 


the solar system is support ed by theoretical studies (e.g., 


1998 ; 


Jacobson et al. 


2 014P and isotope chronometry (e.g., 


- - — 

Wetherill 

1985; 

Koku 

30 & Ida 

. Touboul et al. 

2007; 

Hallidav 


2008). Therefore, these observations indicate that the giant impact stage during which a 


few Earth-sized planets form around 1 AU is not very common around ot 


10%, which is consistent with the estimate made by 


Jackson & Wyatt 


rer stars, probably 


( 2012 ). 


Because many extrasolar planets have been discovered so far, we can compare our 
results with the planet occurrence rates of extrasolar systems. In particular, thanks to the 
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Kepler space telescope, the occurrence rates of Earth-sized planets and super-Earths have 
been recently investigated. Although there is not sufficient data about Earth-sized planets 
and super-Earths located around 1 AU, recent Kepler data show that the occurrence rate 
of planets with radii of 1 R e < R < 2 R ffi and with orbital periods of 150 < P < 250 days 


is 9.6 ± 3.3% for MKGF stars ({Mulders et ah 20 15). where R @ is the radius of the Earth. 


This reported rate seems to be consistent with our speculation that the giant impact stage 
as experienced in our solar system is not very common. 


Because many short-period super-Earths have 
of these planets are reliable. By using Kepler data, 


Deen discoverec 


Petigura et ah 


, the occurrence rates 


(2013) found that the 


occurrence rate of planets with radii of 1 R® < R < 2and with orbital periods of 
5 < P < 100 days is 26 ± 3%. Here, we consider the case of giant impacts that make 
super-Earths located at less than 0.4 AU (which corresponds to an orbital period of 100 days 
around a G-type star). Because the orbital period at 0.4 AU is shorter than that at 1 AU, 
the orbital instability among protoplanets that causes giant impacts takes place very quickly. 
Therefore, the timescale of the giant impact stage for these planetary systems is much 


shorter than that of our terrestrial planets. Moreover, to m ake super - Eart 


protoplanets should be larger than Mars-sized protoplanets flKokubo fe Ida 


is in situ, the 


protoplanets also shorten the timescale of the giant impact stage (K o ku bo et al 


2002 


. Larger 


2006). 


Therefore, even if giant impacts take place in the systems of short-period super-Earths, 
the duration of the giant impact stage should be much shorter (~ 1 Myr) than that for 
our terrestrial planets (~ 100 Myr). Moreover, the lifetime of debris disks produced by 
giant impacts should be also very short, because it strongly depends on the distance of 
the debris disk from the central star (see Eq. ([2])), although the effect of large Mi rg , A a 
and e on prolonging the lifetime of debris disks should be considered. The lifetime of a 
debris disk around 0.4 AU is 50 times shorter than that around 1 AU. Therefore, the 
probability of observing a giant impact stage for short-period super-Earths is extremely low. 
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One additional point that might support non-detectability of super-Earth giant impacts 
is that collisions between super-Earths likely produce more vaporized ejected materials, 
because of their larger escape velocities and higher impact velocity. Because the size of 
condensed materials would be very small, the lifetime of the debris disks produced by 
collisions between super-Earths would be shorter than that produced by collisions between 
Earth-sized objects. Therefore, it is likely that observed warm debris disks are the result of 
the giant impact stage of our solar system-like terrestrial planet formation. 

HG and HK gratefully acknowledge support from Grant-in-Aid for Scientific Research 
(B) (26287101). This work was also supported by Research Grant 2015 of Kurita Water and 
Environment Foundation. Numerical simulations were in part carried out on computers at 
the Center for Computational Astrophysics, National Astronomical Observatory of Japan. 
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Fig. 1.— The total mass of the ejected fragments M e j e by a single giant impact between 
protoplanets of the same size with a mass of 0.1 M® as a function of impact velocity t'i mp 
and angle 6. The color contour represents M e j e normalized by the total mass of the colliding 
protoplanets M sys . From bottom to top white contours represent M e j e /M sys = 10 -3 ’ 5 — 
IO -0 ' 5 with an interval of IO 0,5 . The impact velocity is normalized by the two-body surface 
escape velocity v esc . A head-on collision correspon ds to 6 = 0. Crosses represent the impact 


conditions of SPH calculations performed by 


GendaetaL 


( 20121 ) 
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Run 


Fig. 2.— Total ejected mass normalized by the Earth’s mass M® during the giant 
impact stage for each run. 
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Fig. 3.— Snapshots of a giant impact for the 11th collision in Run 1. This simulation used 
100,000 SPH particles. After the first contact, the protoplanets escape from each other and 
are no longer gravitationally bound. In this collision, a large amount of fragments (0.05M® 
in total) is ejected. The mass of the largest fragment excluding protoplanets is composed of 
approximately 2,500 SPH particles, which corresponds to 1/160M e (half of the lunar mass). 
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Fig. 4.— Total ejected mass M eje and the mass of largest fragment M lrg produced by each 
giant impact occuring in Runs 1, 21, and 26. 
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Fig. 5.— Evolution of the mass of the debris disk for Runs 1, 21, and 26. Each spike 
corresponds to a giant impact, and the decrease in M disk is caused by collisional cascade. 
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Fig. 6.— Calculated observable 24 /jrn flux scaled by the stellar 24 /jm flux as a function 
of time. The observable flux F° hs is the summation of the debris disk’s flux F,? lsk and the 
stellar flux F*. Observational data of infrared excess are from 
labeled with numb ers 1 and 2, 


numbers 4 to 8, 


Me: 


Carpenter et al 


is et al. 


(120101 ) for number 3, 


Chen et al. 


(20111) for s tars 


Zuckerman et ah 


( 20091 ) for numbers 9 and 10, and 


(1201 111 for 


Beichman et al 


(2006) 


for numbers 11 and 12. 



























































